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1. Introduction 

The discovery of a number of double-charmonium production processes at B factories about a 
decade ago [|l|, ^, has triggered the long-lasting interest. Experimentally, this production mech- 
anism offers a unique environment to search for the new C-even charmonium states, particularly 
those X, Y, Z states, by fitting the recoil mass spectrum against the J/y (y'). The famous examples 
are the establishing of the X(3940) and X (41 60) states along this avenue. 

On the theoretical side, the double-charmonium production provides a powerful arsenal to 
strengthen our understanding toward perturbative QCD, especially toward the application of the 
light-cone approach [Q, ^ and the nonrelativistic QCD (NRQCD) factorization approach []6j] to 
hard exclusive reactions involving heavy quarkonium. 

Thanks to the very clean J/y(y') — > signals, experimentally it is most favorable to re- 
construct those double-charmonium events that involve a J/y (y') meson. Thus far, the most 
intensively-studied double charmonium production process is e + e~ — > J jy + T] C . The lowest-order 
(LO) NRQCD predictions to this process [0, |[ was about one order of magnitude smaller than 
the original Belle measurement [jj]]. This disquieting discrepancy has spurred a great amount of 



theoretical investigations in both NRQCD and light-cone approaches QIOQ . One crucial ingredi- 
ent in alleviating the discrepancy between the NRQCD prediction and the data is the positive and 
substantial next-to-leading order (NLO) perturbative corrections [11, 12] l . By contrast, owing to 



some long-standing theoretical obstacles for the helicity -flipped process, e.g. the "endpoint singu- 
larity problem", the 0(a s ) correction to this process has never been worked out in the light-cone 
approach. Therefore, despite some shortcomings, the NRQCD factorization approach remains to 
be the only viable formalism to tackle double-charmonium production which is both based on the 
first principles of QCD and also amenable to the systematical improvement. 

In this talk, we review some recent investigations on the higher order corrections to the double- 
charmonium production processes e + e~ — > J/y+Xc0,ifi( r lc2) at B factories [15, 16, 17, 18, |l9| ]. 



In particular, we address the impact of the 0(a s ) [or 0(a s v 2 ), where v denotes the characteristic 
velocity of the charm quark inside charmonium] corrections on these processes in NRQCD factor- 
ization approach. The theoretical predictions will be confronted with the measurements whenever 
possible. We also remark on a novel theoretical phenomenon, i.e., the occurrence of the double 
logarithms in the NRQCD short-distance coefficients at NLO in a s . 

The rest of the paper is structured as follows. In Sec. Q we recall the helicity selection rule 
relevant for the exclusive double charmonium production processes considered in this work. In 
Sec. [3[ we sketch some key techniques underlying the 0{a s ) calculation. In Sec. fij we consider 
the 0{a s v ) correction to the process e + e~ — > J /y + r\ c , presenting the asymptotic expressions 
of the 0(a s ) NRQCD short-distance coefficients through relative order- v 2 . The phenomenological 
impact of this new correction is also explored. In Sec. |5| we briefly discuss the effects of the 0(a s ) 
corrections to the processes e + e~ — > J '/ 'y ' + ^0,1,2(^2) > and propose that these types of processes 
may be used to unravel the quantum number of the famous X(3872) meson. Finally in Sec. ^, we 
present our personal perspective on those most important problems remaining in this field, which 
urgently await the exploration. 



'The tree-level relativistic corrections to this process have also been investigated ^ |l3[[lj]. Including these cor- 
rections appears to be helpful to further reduce the gap between the NRQCD prediction and the data. 



2 



Some recent progress in understanding exclusive double charmonium production at B factories 



YuJia 




Figure 1: One sample LO diagram and five sample NLO diagrams that contribute to y* -^J/y+H, where 
H is a C-even charmonium state such as r\ c , Xcd,\,2 or r\ C 2. 



2. Helicity selection rule for double-charmonium production 

For the exclusive double-charmonium production processes, it is most informative to look into 
the polarized cross section, where the power-law scaling in each helicity configuration is governed 
by the famous helicity selection rule (HSR) [|2C|]. In the hard-scattering limit 3> m c S> Aqcd {yfs 
stands for the center-of-mass energy of the e + e~ collider, m c for the charm quark mass, and Aqcd 
for the intrinsic QCD scale), the HSR implies that the asymptotic behavior for the rate of producing 
J/\j/ together with a C-even charmonium H (with the leading cc Fock component possessing the 
quantum number 25+1 L) ) in a definite helicity configuration is [0] 

a [e+e-^J/Y(h)+Hfa)] 6+2L (*t \ 2+]h+ ^ 
a[e+e-^n+Ll-} W / ' 



where X\, X2 represent the helicities carried by the J/\j/ and H, respectively Eq. (IA) implies 
that the helicity state which exhibits the slowest asymptotic decrease, thus constitutes the "leading- 
twist" contribution, i.e., o ~ \/s\ is (Ai,A 2 ) = (0,0). For the processes exemplified by e + e — > 
J/W+ f]c{%c\,T}c2), parity invariance forbids the occurrence of the (0,0) configuration, therefore 
they are entirely of the "higher twist" (helicity-flipped) nature. 

3. Techniques in calculating the 0(a s ) corrections 

The NRQCD short-distance coefficient can be most easily obtained by computing the on-shell 
quark amplitude 7* — > cc(fi, 3 Sp) + cc(P2, 2S+l Ly ), with the aid of the covariant spin projec- 
tors [|2l], 0]. At LO in a s , there are 4 diagrams for this parton process; while at NLO in a s , there 
are 20 two-point, 20 three-point, 18 four-point, and 6 five-point one-loop diagrams. Some typical 
diagrams are shown in Fig. [[} 

In spite of resorting to perturbative matching method, a much more economic way is to directly 



extract the short-distance coefficients following the philosophy of threshold expansion [|22J]. That 
is, when projecting each cc pair onto the desired orbital-angular-momentum state, one first expands 
the amplitude in powers of the relative quark momentum q prior to performing the loop integration. 

Consequently, a technical complication arises, that one inevitably encounters some unusual 
one-loop integrals that contain the propagators of (up to) cubic power, due to taking the derivative 
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over q. The Mathematica packages FIRE p3]] and the code Apart $M are utilized to reduce 



these unconventional higher-point one-loop tensor integrals into a minimal set of masters integrals. 
Thanks to the integration-by-part algorithm and partial fraction technique built into these codes, it 
turns out that all the encountered master integrals are nothing but the ordinary 2-point and 3-point 



one-loop scalar integrals, whose analytic expressions can be found in Appendix of Ref. []12[]. 

When adding the contributions of all the diagrams, and after renormalizing the charm quark 
mass and the QCD coupling constant, we end up with both UV and IR finite NLO expressions for 
the 0(a s ) amplitude associated with y* — > J /y + Xc0,\,2{i]c2)- hi comparison, the amplitude for 
7* — > J/Y+ f] c contains an uncanceled IR pole at 0(a s v 2 ). Fortunately, it can be factored into the 
relative order-v 2 NRQCD matrix element via the pull-up mechanism. In either case, one ends up 
with the IR-finite 0(a s ) NRQCD short-distance coefficients for all the helicity amplitudes affiliated 
with e + e~ — > J/y+H. Our calculation explicitly confirms the assertion made in Ref. jj25l], that the 
NRQCD factorization holds beyond tree level for the exclusive production of a 5-wave quarkonium 
plus any higher orbital-angular-momentum quarkonium in e + e~ annihilation. 

4. 0(a s v 2 ) correction to e + e~ J/\j/ + r) c 

The production rate for e + e~ — > J/y+ Tf c can be expressed as 

4na 2 /|P|\ 3 ,„, m2 



a[e + e- ^ J/ v+r]c ] = — \l-±) \ G (s)\\ (4.1) 

where |P| signifies the magnitude of the momentum carried by the J /iff (rj c ) in the center-of-mass 
frame. G(s) is the 7/y/ + rj c time-like electromagnetic form factor, defined through (J/\ff(Pi, A) + 
77 c (P 2 )|./emI > = iG{s)e^ pa Pi v P 2p el{X), where 7^ M is the electromagnetic current, Pi (£ ff (A)) 
denote the momentum (polarization vector) of the J /y, and P2 the momentum of the 7] c , respec- 
tively. 

4.1 NRQCD factorization formula 

NRQCD factorization allows one to factorize the J /\y + r] c electromagnetic form factor as 



G(s) = J4Mj /xif M Vc {J/\ir\^<r-ex\0}{T] c \Y^x\0) [co + c 2 ,i(v 2 ) 7/v/ + c 2 ,2(v 2 ) t?c + ---] ,(4.2) 



where Co and C2j are the dimensionless short-distance coefficients that depend on m 2 /s. For sim- 
plicity, we have also introduced the following dimensionless ratios of NRQCD matrix elements 
to signify the <9(v 2 ) corrections: (v 2 ) J/¥ = (J/\jf(l)\Y^(-^ D) 2 <T ■ e(A)^|0)/(m 2 {J/y(X)\y^O- 

B(X)X\0)), (v 2 ) Vc = ^cWK-^) 2 Xm{m 2 c ^ c \^x\0)), where ^ T) X = W^X ~ {^wV X- 
Inserting Eq. into (|Q1), one can decompose the cross section into the 0(v°) and 0(v 2 ) 



pieces, o[e + e ->• J/y+ T] c ] = a + o 2 + O(ov 4 ), where 

ct = 8 " a2 " , - ( 3'""' >3/2 (^) J / t ,(^),.ko| 2 , (4.3a) 
a, = * a ^ ( 3 '-* ) V w*>» (43b) 
I (n"^-^"!" + 4Re l CoC 5.']) <" 2 >J/v 4" ( : ^7kor' + 4Re[co42]')(v 2 ),, r L 
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We have employed the Gremm-Kapustin relation [^] to eliminate the explicit occurrences of Mjiy 



and Mrj c in (43). To condense the notation, we have introduced the following symbols: r = 4ml /s, 



m J/v = \(J/w(e)\¥ l <y-ex\0)\ 2 , and (0^ = |(t? c |^|0)| 2 . 

4.2 Various NRQCD short-distance coefficients 

We organize the coefficients c, (i = 0,2) in power series of the strong coupling constant, i.e., 
Ci = cf^ + %c- + • • • . Accordingly, one may decompose the cross section a, into o- + G-^ 
(i = 0, 2) as well. Thus far, the only missing piece is . The tree-level short-distance coefficients 
through 0(y 2 ) have been available long ago [0]: 

(0) _ 327^^ jo) _ 3-10r (Q) (0) _ 2-5r (0) 

~ N c m c s 2 ' 2' 1 " 6 ' 2 < 2 " 3 ' { ' 

7 — 1 4 

where e c = | is the electric charge of the charm quark, and Cp = = |. 

The 0{a s ) NRQCD short-distance coefficients and c^ 1 ,- are generally complex- valued, 
and cumbersomely-looking. Nevertheless, it is much more illuminating to look at their asymptotic 
expressions in the limit y/s 3> ra c : 

J'){^r\ _m J r J 1, ' 5\ , /13,, 5,.,,.. 41 



\ / asym v 

31 , 1 , 2 „ 893, „ 5 2 497\ (\ „ 3, 19 



+ _ lnr __ W2 + _ ln2 __^__ .+„ _ ft + _ lnr+ _ ln2+ ^ 1 J, (4 . 5b) 





asym 



1 11 2 241 1 2 99 \ /l „ 1 11 1 v I 

■mr-— ln 2 2 + — In2--7T 2 -— )+in \ -ft + rlnr+ — ln2 - — (4.5c) 



24 8 144 8 16/ \ 4 6 12 24 

where j8o = yCi — |n/ is the one-loop coefficient of the QCD /3 function, and ft/ = 4 denotes the 
number of active quark flavors. /x r denotes the renormalization scale, and \lf signifies the NRQCD 
factorization scale in the MS scheme, which naturally ranges from m c v to m c . 



2 (1) 

As first pointed out in Ref. [27], a peculiar double-logarithmic correction o< hi z r arises in c 



for this helicity-fiipped process, and our (4.5a) exactly agrees with the corresponding expression 
there 2 . Eqs. ( |4.5[ ) imply that the double logarithms survive at 0(a s v 2 ) as well. With y/s S> m c , 
one presumably needs to sum these types of logarithms to all orders in a s to warrant a reliable 



prediction. Such a resummation may even be mandatoiy at the fi-factory energy Q27|]. At present, 
how to fulfill this goal remains to be a thorny challenge. 



2 Our result slightly differs from Ref. [jl^| on the imaginary part of Cq , though it does not affect the phenomenology 
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Table 1: The individual contribution to the predicted a[e + e — > J/y/+ rj c ] at yfs = 10.58 GeV, as specified 
by the powers of a s and v 2 . We take a(s/s) = 1 /130.9, = 0.338 GeV, m c = 1 .4 GeV, and \i f = m c . The 
LO NRQCD matrix elements are ~ {@\)r\ c = 0.387 GeV 3 . The Gremm-Kapustin relation is used 

to obtain (v 2 ),//^ = 0.223 and (v 2 )^ = 0.133. The cross sections are in units of fb. 
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Figure 2: The fi r - and -^-dependence of the cross section for e + e — > J/y+ T\ c - The 5 curves from bottom 
to top are (solid line), <7q ' 



line), and <x 



(0) 



,0) 



?P + a (1) (dashed 

ai 1 ^ + ai 1 ' (solid line), respectively. In the left panel, we fix = 10.58 GeV, and 



(i) 



a) (dashed line), (7^°' + CTq 1 ' (solid line), 0^ 



the blue and green bands represent the measured cross sections by the Belle and BaBar experiments, with 
respective systematic and statistical errors added in quadrature. 



4.3 Phenomenological impact of the 0((X s v 2 ) correction 

Apart from the ambiguity in the values of m c and NRQCD matrix elements, the freedom 
of choosing the scale entering the strong coupling constant leads to a large uncertainty for our 
prediction to o[e + e~ — > J f] c ]- This is a serious drawback of the NRQCD approach [j28|]. For 
simplicity we assign all the occurring a s with a common scale, and choose \i r = \fs/2 and 
H r = 2m c , respectively, hoping that the less biased results interpolate somewhere in between. 

Table [j] lists the predicted o[e + e~ — > J/y + i} c ] with two sets of \l r , organized in double 
expansions of a s and v. We reproduce the well4mown results, i.e., the positive and substantial 
0(a s ) correction [ pT| , ^2Q , and the positive but less pronounced 0(v 2 ) correction [0, |l3, 14]. Taking 
r = 0.07, which is relevant for the B factory energy, Refc^ 1 ,- /c^ } ] (i = 1 , 2) turn out to be both large 
and negative. One might naively expect that including the new 0(a s v 2 ) correction would largely 
dilute the existing 0(y 2 ) term. However, the new correction to the cross section, a!} , is actually 
positive and modest. This may be attributed to the accidental cancelation between the two terms in 
the prefactor of (v 2 )h in ( 4.3b ), which represent two difference sources of relativistic correction. 

In Fig. ||, we plot the o[e + e~ — > J/y+Vc] as a function of \i r and yfs. At the B factory energy, 
incorporating the new correction appears not to make a big difference. When \L r is relatively 
small, the state-of-the-art NRQCD prediction converges to the BaBar measurement within errors. 
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Had we taken somewhat larger values of the NRQCD matrix elements {&\)u as in [13, Qj, the 
agreement with the B factory measurements would improve. 

As shown in the right panel of Fig. ||, o[e + e~ — > //y + He] drops steeply as \fs increases, 
reflecting the helicity-flipped nature of this process, a ~ \/s 4 . In contrast to its minor impact at B 
factory, the 0(a s v 2 ) correction turns out to be much more relevant at higher -y/s. 

5. 0(a s ) correction to e + e = -►.//> + 1,2(^2 ) 

The NLO perturbative corrections to the double-charmonium production processes e + e~ — > 
J / Y + Xc0,i,2{ilc2) have also been investigated recently [|TJ, |TJ, |l^]. The calculational technique 
resembles that for the 0(a s v 2 ) correction to e + e~ — > J/y+ r\ c , and has been reviewed in Sec. ||. 
The processes e + e~ — > J /\ff + ^-0,1,2 have been independently studied by two groups, and they 
agree with each other [jTJ], |T|] . 

The 0(a s ) correction to e + e~ — > J/y + XcO appears to be substantial. With some reasonable 
input parameters, and the renormalization scale \x r set as including the 0(a s ) correction en- 
hances the LO cross section from 4.8 fb to 8.6 fb. This value is compatible with the Belle and 
BaBar measurements within errors. However, the NLO perturbative prediction to the + XcO 



production rate is still significantly below the central value of the Belle measurement [15, |lq ]. 

The situation for e + e~ — > J j\\f + Xd,2 is less clear. The impact of the NLO perturbative cor- 
rections to these processes seems to be modest, even with the sign uncertain, depending on the 
different choices of jx r . The predicted cross sections for both processes are around 1 fb, which 
are almost one order of magnitude smaller than that for e + e~ — > J/y + XcO- Encouragingly, the 
recent Belle experiment [ p9| ] did observe a considerable number of e + e~ — > J/y + Xd,2 events, 
from which one may roughly estimate the corresponding production cross sections. They appear 
to be qualitatively consistent with our expectations. 

Recently there has arisen some controversy about the canonical charmonium option of the 
X(3872) meson, whether it being T] C 2 or %' cl [p0|]. Motivated by this concern, we have also per- 
formed a comparative study for e + e~ — > J r\ C 2 an d e + e~ — > J /y + x' c \ at B-factory energy, 
hoping that it may provide some guidance to unravel the quantum number of the X (3872) meson in 
the future experiments Jl9|]. The NLO perturbative correction to the former process is of medium 
size. With ji. r taken as implementing the 0(a s ) correction enhances the LO cross section 

from 0.22 fb to 0.29 fb. The production rate of the latter process is about 6-7 times greater than 
the former, thereby it seems realistic to observe the J/y + x' c i signals based on the current 1 ab _1 
Belle data sample, if the x' c \ is indeed the narrow X(3872) meson. 



Refs. [|16J, |19J also conduct a comprehensive study on the polarized cross sections for e + e~ — > 
J/Y+XcO 1 2{f]c2)- It is found that the bulk of the total cross section comes from the (0, ±1) helicity 
channels for e + e~ J/y + Xci, from the (0,0) and (±1,0) helicity states for e + e~ J/y + Xc2, 
and from the (±1,0) states for e + e~ — > J/y+rj C 2- The hierarchy among the various helicity 
channels appears to often conflict with what is expected from the HSR. It will be interesting for the 
future experiments to concretely test these polarization patterns. 

By working out the asymptotic expressions of the various helicity amplitudes for the processes 



7* — > J ' 1 'y/ '+ XcQ,\,2{f]c2), we firmly confirm the pattern speculated in Ref. [27]: The hard exclusive 



reaction involving double charmonium at leading twist can only host the single collinear logarithm 



7 



Some recent progress in understanding exclusive double charmonium production at B factories 



YuJia 



Ins/ml at NLO in a s , while the double logarithms of form In 2 s/ra 2 are always associated with 
those helicity-suppressed channels. 

6. Outlook 

After a decade of intensive study, our understanding of exclusive double charmonium produc- 
tion has gradually matured. The most notable lesson is perhaps that, the NRQCD factorization 
approach has proved to be a successful and indispensable tool in dealing with hard exclusive re- 
actions involving heavy quarkonium. However, in our opinion, this research area is still far from 
being closed, and there remain some important questions to be answered. In the following, we 
enumerate two topics which may urgently beg for the exploration. 

We have reviewed some recent advances in the NLO perturbative correction to double char- 
monium production processes. Aside from e + e~ — > J /\y + r\ c , the relativistic correction has hardly 
been investigated for any double-charmonium production process involving the P, D-wave charmo- 
nium. There is no ground to believe their effects are less important than the perturbative corrections. 
A practical difficulty to assess the relativistic correction for these processes is that, in general more 
NRQCD matrix elements than for the 5-wave charmonium will come into play, about whose val- 
ues we have almost no any clue. Hopefully, the lattice NRQCD simulation will eventually provide 
some useful information for those long-distance matrix elements. 

A great theoretical challenge is to tame the double logarithms of form In 2 {s/m 2 c ) in the 0(a s ) 
NRQCD short-distance coefficients which are always affiliated with the helicity-suppressed ex- 
clusive quarkonium production channels. The occurrence of these process-dependent, (positive) 
double logarithms severely jeopardizes the reliability of the fix-order perturbation theory predic- 
tion. We note that some important progress has been made recently in tracing the origin of these 
double logarithms at one-loop order (differentiating the harmless Sudakov double logarithm from 
the problematic endpoint double logarithm) [pl|]. Nevertheless, there is still a long way to go to 
finally develop a systematic control over these endpoint double logarithms, e.g., to resum them to 
all orders in a s . In our perspective, the occurrence of these double logarithms is likely intertwined 
with the long-standing failure in applying the light-cone approach to the hard exclusive reactions 
beyond tree level. Looking on the bright side, with explicit expressions of the 0{a s ) NRQCD 
short-distance coefficients for many channels at our disposal, one may view the exclusive double 
quarkonium production as a fertile theoretical laboratory, from which some fresh insight may be 
gained by reexamining those old problems of the light-cone approach. 
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